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ABSTRACT 

Context. In many plasmas, long-lived metastable atomic levels are depopulated by collisions (quenched) before they decay radiatively. 
In low-density regions, however, the low collision rate may allow depopulation by electric dipole (El) forbidden radiative transitions, 
so-called forbidden lines (mainly Ml and E2 transitions) . If the atomic transition data is known, these lines are indicators of physical 
plasma conditions and used for abundance determination. 

Aims. Transition rates can be derived by combining relative intensities between the decay channels, so-called branching fractions 
(BFs), and the radiative lifetime of the common upper level. We use this approach for forbidden [Sc n] lines, along with new calcula- 
tions. 

Methods. Neither BFs for forbidden lines, nor lifetimes of metastable levels, are easily measured in a laboratory. Therefore, astrophys- 
ical BFs measured in Space Telescope Imaging Spectrograph (STIS) spectra of the strontium filament of Eta Carinae are combined 
with lifetime measurements using a laser probing technique on a stored ion-beam (CRYRING facility, MSL, Stockholm). These quan- 
tities are used to derive the absolute transition rates (A-values). New theoretical transition rates and lifetimes are calulated using the 
CIV3 code. 

Results. We report experimental lifetimes of the Sen levels 3d 2 a 3 Po,i,2 with lifetimes 1.28, 1.42, and 1.24 s, respectively, and tran- 
sition rates for lines from these levels down to 3d4s a 3 D in the region 8270-8390A. These are the most important forbidden [Sell] 
transitions. New calculations for lines and metastable lifetimes are also presented, and are in good agreement with the experimental 
data. 

Key words. Atomic data - Radiation mechanisms: general - Methods: laboratory - Stars: individual: Eta Carinae 



. — 1- Introduction The FERRUM project was initiated to obtain experimental 
, data and evaluate theoretical data for the iron group ele ments 
'jjj . Physical conditions and elemental abundances of stellar and neb- used in astrophysical applications dJohansson et al.ll2002l) . This 
03 | ular objects are derived from quantitative analysis of spectra. For i s an international collaboration with the aim to produce reli- 
H II regions and gaseous nebulae, emission lines are used to de- able oscillator strengths and transition rates for transitions used 
rive these properties based on simple line ratios or more detailed i n astrophysical analysis. One branch within this collaboration is 
modeling. In either case, accurate atomic data is needed to obtain to measure metastable lifetimes and, in combination with astro- 
reliable determinations. One important property is the radiative physical branching fractions, derive absolute transition rates for 
transition rate for spontaneous decay. forbidden lines. Forbidden lines are, by definition, fines that do 
For complex spectra, which have many lines, not all transi- not occur by electric dipole radiation (El ), but rather with higher 
tion rates can be measured. The majority of the transition data multipole radiation, in particular magnetic dipole (Ml) or elec- 
need to be determined through calculations. Experimental data trie quadropole (E2). Parity forbidden lines are transitions be- 
with reliable uncertainties are measured for important transitions tween levels of equal parity, and occur between the lowest levels, 
and used to compare with the calculations. which cannot decay via El transitions. Forbidden lines with even 

higher multipole radiation can have transitions between levels of 
opposite parity. In this case, it is often the change in /-value that 
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* Based on observations made with the NAS A/ES A Hubble Space radiative lifetimes of the order of seconds, compared to nanosec- 

Telescope, obtained at the Space Telescope Science Institute. STScI is onds for normal levels, 
operated by the Association of Universities for Research in Astronomy, 

Inc. under NASA contract NAS 5-26555. Laboratory data is obtained Due to the long lifetimes of metastable levels, collisions are 

at the CRYRING facility at the Manne Siegbahn Laboratory (MSL), often the dominating deexcitation process in dense plasmas. In 

Stockholm University, Sweden. low-density plasmas, as in gaseous nebulae and HII regions, the 
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Fig. 1. Partial energy level diagram of Sc n, showing the three 
lowest configurations. The lines measured in the present paper 
are 3d4s a 3 D - 3d 2 a 3 P. 



radiative decay can be of similar strength or the main process. 
This sensitivity to density makes the forbidden lines good for 
density diagnostics, and also as temperature probes by simple 
line ratio measurements. The same collision effect occurs in the 
laboratory, and a low pressure environment is required to mea- 
sure the lifetimes of metastable levels. In the present work, we 
used a storage ring with ultra-high vacuum where the ions can be 
stored for a long enough time to allow for measurement. Another 
consequence of the sensitivity to collisions is the problem of pro- 
ducing forbidden lines to measure branching fractions, i.e. the 
fraction of decay in the different decay channels from a specific 
upper level. We use HS T/STIS spectra of Eta Carinae to mea- 
sure the lines. 

Eta Carinae is one of the most massive stars known in the 
Galaxy. It has a violent history and has experienced several 
nonterminal explosions of which the so-called Great Eruption 
in the 1840's produced its today well- known bipolar lobes and 
the intervening disklike structure (e.g. [Davidson & Humphreys! 
1 19971; [Smith et alj20 04). One region of the disk shows a remark- 
able spectrum, among others forbidden [Srn] along with many 
permitted and forbidden lines of iron group elements, specifi- 
cally Ti ii ( iHartman et al. 2004). This region was given the name 
'Strontium Filament' (SrF) from the identification of [Srn], 
which previously was not observed in any othe r astronomical 
objec t (Zethson et al. 2001). In a previous paper Hartma n et alj 
(2005) derived A-values for [Tin] with these spectra. Photo- 
ionization modelling indicate peculiar abund ance pattern s and 
photo-excitation in this region (Bautista et al. 2002, 2006). 

In the spectrum of the SrF of Eta Carinae (Hartman et al.l 

12004 forbidden l ines of f Sc n] w ere observed, 3d4sa 3 D-3d 2 a 3 P 
(Multiplet 3F of lMoord dl945l) ). We measured the branching 
fractions for these lines to determine the transition rates. In 
Fig.Q] a partial energy level diagram of Se n is shown, based on 
the analysis by iJohansson & Litzenl JT980). Included are levels 
from the three lowest configurations: 3d4s, 3d 2 , and 4s 2 . Lines 
between these levels are parity forbidden. 

In Sec. [2] we describe the experimental technique: the life- 
time measurements and the derivation of astrophysical branch- 
ing fractions. In Sec. [3] we present new calculations using the 
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Fig. 2. Energy level diagram showing the pumping channels and 
major fluorescence decay used in the present experiment, with 
their approximate wavelengths. The laser is tuned to match the 
transitions to z 3 D or z 3 P, and the subsequent decay to a 3 D is 
observed. The level energies are not to scale. 

CIV3 code and in Sec. |4] the results are discussed, along with 
uncertainties. 



2. Experimental technique 

A widely-used technique to derive absolute transition rates (A), 
or oscillator strengths, is the combination of the lifetime (r) of 
an upper level and the branching fractions (BFs) of the different 
decay channels from the same level: 



A,d = 



where 



BF„i = 



BF„ 



1,1 A* 



and 



T u = 



(1) 



(2) 



(3) 



The subscripts u and I represent the upper and lower level, re- 
spectively. The technique is primarily used for electric dipole 
transitions (El). In that case, the lifetime is often measured with 
a laser induced fluorescence (LIF) technique, and the BFs are 
measured in a laboratory light source. For forbidden lines, this 
technique cannot be used since the timescale is totally different: 
seconds, instead of nanoseconds. 

2. 1 . Measurement of metastable lifetimes in Sc n 

At the CRYRING facility, a technique has been developed to 
measure lifetimes of metastable levels by laser probing of a 
stored ion beam. The ultra-high vacuum (~10 -11 mbar) makes 
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Table 1. Lifetimes in seconds. 
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Fig. 3. Decay curve for the a P2 level at base pressure. 
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* The uncertainty quoted is a one standard deviation of the linear fit to 
the Stern- Vollmer plot. 
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Fig. 5. Comparisons of experimental and theoretical lifetimes for 
the 3d 2 a 3 P levels from the present wor k, and to previous calcu- 
lations bv lWarner & K irkpatrickl (|l969l) . 
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Fig. 4. Stern- Vollmer plot for the a 3 P2 level. The apparent life- 
time is measured at base pressure and raised pressure, and the 
true radiative lifetime is derived by extrapolation to zero pres- 
sure. 



collisons rare, and the metastable ions survive for a time compa- 
rable to several lifetimes. The laser probing technique (LPT) and 
its later developments have been decribed in a number of publi- 
cation s dLidberg etalJI 19991 lMannervikll2003t iMannervik et alj 
120051) . Here we only discuss the general idea of the technique 
and issues specific to the present measurement, and refer to pre- 
vious publications for a more detailed description. 

The Sc + ions are produced in an ion source and, after be- 
ing accelerated to an energy of 40 keV and passing through an 
isotope-selecting bending magnet, are injected into the storage 
ring. A continuous wave (cw) tunable dye laser, with a mechan- 
ical chopper, is used to probe the stored ions. The wavelength 
is tuned to match a strong transition from the metastable level 
to an upper level with opposite parity (Fig. |2). The subsequent 
decay is fast, a few nanoseconds. The intensity of the prompt 
decay, the fluorescence, is proportional to the number of ions 
in the metastable level probed at a certain delay after injection 
into the ring. By varying the delay time of the laser pulse after 
injection into the ring, the time dependence of the population in 
the metastable level is monitored. For each injection, the popula- 
tion at one time is measured and corresponds to one point on the 



decay curve. Subsequent injections into the ring are probed at 
different time delays, and in this way the decay curve is built up. 
The technique is destructive in the sense that once a laser pulse 
is applied, the level is emptied and a new set of ions must to be 
injected for the next measurement. It is thus important that the 
number of ions in the level measured do not vary too much and 
that the variations can be monitored. During the measurement, a 
number of properties, such as the beam lifetime, number of ions 
in each injection, and the fraction of ions in the level of inter- 
est, are monitored to allow for c orrections of pos sible drifts and 
changes in the ion production (Mannervik 2003])- A corrected 
decay curve is shown in Fig. [3] 

Despite the ultra-high vacuum, the stored ions are affected by 
collisions with rest gas. Both collisional deexcitation (quench- 
ing) and excitation (repopulation) are observed and corrected 
for. A detailed description of the corrections is given in previous 
publications (see above). The quenching is observed as a pres- 
sure dependence in the measured lifetime. The collisions add a 
decay rate in addition to the radiative rate of interest. This col- 
lisional decay rate can be subtracted by measuring the lifetime 
as a function of pre ssure and ex trapolating to zero pressure in a 
Stern- Vollmer plot dDemtroderll 1 9961) . see Fig. [4] 

The experimental lifetimes of 3d 2 a 3 Po,i,2 are presented in 
Table 1 along with uncertainties, which are one sigma uncer- 
tainties to the linear fit of the Stern- Vollmer plot. 



2.2. Branching fractions for3d4s a 3 D-3d 2 a 3 P 

With the high spatial resolution of the STIS spectrograph on- 
board HS T, the different parts of the nebulosity and the ejecta 
around Eta Carinae can be spatially resolved. Many of these 
have low densities, and allow for parity forbidden lines to be 
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Fig. 6. Spectrum showing the [Sen] lines used in the present 
analysis. The upper level of each group is indicated. The tran- 
sitions are listed in Table [3] 

produced. One region rich in forbidden lines is the so-called 
Strontium Filament. 

Among the many forbidden iron-group elements lines, the 
a 3 D-a 3 P transitions in [Sen] are observed, as shown in Fig. [6] 
Other forbidden lines of [Sen] fall longward of the observed 
wavelength region. Lines from a'G^ with energy similar to the 
a 3 P levels, are masked in the spectrum by other strong lines. 
The major decay from b'So would have a wavelength around 
4270 A, but this line is not observed, probably due to high exci- 
tation energy and thus low population of the b'So level. A prob- 
lem that occurs when using astrophysical spectra instead of lab- 
oratory sources to measure BFs is the blending of lines from 
other elements. The present lines fall in the region where the HI 
Paschen series has numerous lines, but from other spatial regions 
of the nebula. However, theser are weak and can be avoided. The 
8263 A [Sc n] line is corrected for blending with an [Fe n] line, 
which contributes to the observed emission feature. The amount 
of blending is determined to be 10% of the total intensity, and is 
derived using other [Fe n] lines, that are strong in other spatial 
regions along the line-of-sight. The a 3 Pi 2 each have two major 
transitions each, whereas a 3 Po has one single dominant decay 
channel. Astrophysical spectra need to be corrected for interstel- 
lar reddening, but since the lines span a very limited wavelength 
region, 8250-8380A, this is not necessary for the present lines. 

3. Calculations 

We have undertaken an extensive calculation of rates of forbid- 
den transitions in Sc n using configuration interaction wave func- 
tions, which are expressed in the form 

M 

V(LSJ) = Yj a &( a ' LSJ ) W 

where {cD, } are single-configuration functions (CSFs) and the ex- 
pansions, in general, include summations over L and S, thus 
allowing mixing between different LS terms with the same J 
value. For a specific choice of {<]>,}, the expansion coefficients 
{a,} are the eigenvector components of the Hamiltonian matrix 
elements H,j = < <&j\H\<&j >. The relativistic effects that make 
these forbidden transitions possible were incorporated into the 
calculation using the Breit-Pauli approximation, which should 
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Table 2. Configurations used in full calculation 
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be adequate for moderately heavy ions such as Sc n. Specifically, 
we included in our Hamiltonian the one-body mass-correction 
and Darwin terms and the nuclear spin-orbit interaction, as well 
as the two-body spin-other-orbit and spin-spin terms. 

The states involved in this calculation are normally labelled 
3d 2 , 3d4s, or 4s 2 . However, the optimal functions for these states 
differ quite substantially, so, for example, the optimal 3d orbital 
for 3d 2 states is quite different from that for 3d4s states. Our 
calculations assume that we are using a common set of orthogo- 
nal orbitals to describe all the states and we need to allow for this 
effective non-orthogonality by introducing additional, correcting 
orbitals. 

The calculations d escribed here have been undertaken 
with the code CIV 3 (iHibbed I1975L iGlass & Hibbertl 11978b 
iHibbert et aTlll991h . For the ls,2s,2p,3s,3p, and 4s orbitals we 
used the Hartree-Fock radial functions given by Clementi and 
Roetti (1974) for the [ls 2 2s 2 2p 6 3s 2 3p 6 ]3d4s 3 D state. We opti- 
mised 3d on the 3d 2 3 P state rather than 3d4s 3 D to provide a 
balance between different states and because the correcting or- 
bital 4d can be then treated as a linear correction rather than 
quadratic. Finally, 4p,4f, and 5d are correlation orbitals for the 
outer subshell. The full process of optimising the radial func- 
tions and descriptions of the configuration sets will be given in a 
separate paper (Lundin et al., in prep) 

The full set of configurations used, for all LS symmetries, is 
shown in Table [2] They comprise all the configurations that can 
be constructed for valence shell correlation, together with the 
important 3p 2 — >3d 2 correlation effect in the dominant configu- 
rations (3d 2 and 3d4s, together with the effects of the correction 
orbitals that give rise to 3d4d and 3d5d). We have also allowed 
the main effect of the replacement of 3p 2 by two other orbital 
functions. The even parity LS symmetries in our work are 31 D, 
3 P, 'S, 3 F and 'G, appropriately mixed for each /-value. In this 
way, we are able to determine wave functions for each of the 14 
levels normally labelled with 3d 2 , 3d4s, or 4s 2 . 

To improve the calculations, we make small corrections to 
the diagonal elements of the Hamiltonian matrix (Hibbert 1996) 
so that the final cal culated energy separation s agree with the ex- 
perimental values dJohansson & Litz6nlll980l) . We have termed 
this process 'fine-tuning'. We have found that it improves the 
accuracy of the mixing coefficients (a,- in (0]i) and hence the cal- 
culated transition rates. We find that there is very little mixing 
between states of different LS , especially for the triplets studied 
in this paper. 

4. Results 

The measured and calculated lifetimes are given and compared 
in Tableland Fig. [5] Our experimental values are 1 .28, 1 .42, and 
1.24 s for the 3d 2 a 3 P levels, for 7=0,1,2, respectively, and our 
calculated values are 1.44, 1.43, and 1.43 s for the same levels. 
E2 transitions dominate the calcul ation of the lifetimes of the 3d 2 
levels. The previous calculation by War ner & Kirkpatrickl(ll969l) 
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Table 3. The strongest lines from the a 3 P levels. The lines lacking experimental branching fractions (BFs) are either too weak to be 
detected or have wavelengths outside the observed spectral region. The residual is the total contribution from the transitions from 
the upper level not included in the table. 



upper level 


lower level 


X vac (A) 


BF„ P (%) 


BF„ ; (%) 


A eX p (s ) 


Atheory (S ) 


a J P 


a 3 F 2 


13752.83 




4.34 




0.0302 




a 3 D 2 


8328.91 


95.6+4 


95.6 


0.75+0.08 


0.666 




Weaker transitions (total) 


0.06 








a 3 P! 


a 3 F 2 


13701.20 




1.48 




0.0103 




a 3 F 3 


13854.39 




2.79 




0.0195 




a 3 D 2 


8309.94 




7.98 




0.0557 




a 3 D 3 


8386.63 


40.0+6 


43.0 


0.28+0.06 


0.300 




a 3 D, 


8263.44 


47.7+7 


44.7 


0.34+0.07 


0.312 




Weaker transitions (total) 


0.05 








a 3 P 2 


a 3 F 3 


13753.55 




0.87 




0.0061 




a 3 F 4 


13953.56 




3.10 




0.0217 




a 3 D 2 


8273.56 


27.8+5 


34.2 


0.23+0.04 


0.239 




a 3 D 3 


8349.57 


58.7+8 


52.3 


0.47+0.07 


0.366 




a 3 Di 


8227.46 




9.09 




0.0636 




Weaker transitions (total) 


0.04 









"The uncertainty is derived from the uncertainty of the lifetimes and the branching fractions. 



gives 1.07 s for all levels. The Warner and Kirkpatrick calcula- 
tion is based on a single-configuration approximation. This ne- 
glects configuration interaction, which for Sc n can be important. 
The energy levels are obtained semi-empirically (effectively ex- 
perimental transition energies are used) and the mixing between 
different LS terms is therefore approximately correct. The radial 
integrals of the dipole matrix elements are obtained using scaled 
Thomas-Fermi orbitals, and are therefore obtained without the 
use of the variational principle. Hence their results, while being 
of reasonable accuracy, are not so reliable as the present calcu- 
lations. 

The uncertainty for the experimental lifetimes are 4%-15%, 
and presented in Table 1 . The uncertainty quoted is one standard 
deviation of the linear fit to the Stern- Vollmer plot. The vari- 
ations in uncertainty between the three investigated levels are 
mainly due to variations in experimental conditions. We mea- 
sured the lifetime of a 3 P2 separately from the lifetimes of the 
a 3 Po ; i levels. During the study of the a 3 P2 level, more favourable 
experimental conditions were achieved. This included a higher 
ion current in the storage ring as well as higher power from 
the ring dye laser which in turn resulted in better statistics with 
smaller uncertainties. 

The smaller uncertainty associated with the lifetime of the 
3 P2 level is also a consequence of the fact that the we mea- 
sured the lifetime of this level in a wider relative pressure range 
compared to the measurements of the 3 Po,i levels, which makes 
the extrapolation to zero pressure more certain. It should, how- 
ever, be emphasized that in all three measurements, the effect 
of collisional quenching is small since the differences between 
the deduced pure radiative lifetimes and the longest individual 
lifetimes measured at base pressure is 2-1 1% . 

The measured branching fractions are presented with error 
bars in Fig. [7] along with the theoretical values. The a 3 Po value 
is not included since only one line carries the major intensity. 
This means that no experimental BF can be determined for lines 
from this level. The uncertainties of the BFs are in the range 4- 
16% and include contributions from intensity measurement, line 
deblending and correction for unobserved lines (residual). The 
residual is the total BF for all unmeasureable lines: 

i) either falling outside the observed wavelength region, or 



ii) too weak to be measured. 
The residuals for a 3 Po,i,2 are, including the transitions to a 3 F, 
0.05, 0.12, 0.14, respectively, and an uncertainty of 100% is as- 
signed for these values. The agreement with the calculated val- 
ues is within the 1-sigma error bars for most transitions. 

The resulting experimental A-values, derived from the mea- 
sured lifetimes and astrophysical BFs, are presented in Table [3] 
and Fig.|8]and are in the range 0.23 to 0.75 s . The uncertainties 
of the experimental transition rates are in the range 11-21%, and 
result from uncertainties in the lifetimes and in the BFs, includ- 
ing uncertainty of unmeasureable lines (the residual, discussed 
above). The contributions to the total uncertainty from the life- 
time and branching fractions are of similar importance. There is 
relatively little configuration mixing between the states, so we 
would expect the level of accuracy of the calculated values to 
be comparable with that encountered in El transitions for which 
there is little mixing. For the theoretical values, an accuracy of 
5-10% could reasonably be assigned. 

5. Conclusions 

We have measured lifetimes of the 3d 2 3 Po,i,2 levels in Sen. 
By combining these with astrophysical branching fractions from 
HST/STIS spectra form Eta Carinae, experimental transition 
rates are derived for the parity forbidden [Sen] lines 3d4s a 3 D - 
3d 2 a 3 P. The values are compared with new CIV3 calculations, 
and the agreement is good. 
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